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Abstract 

Investigation of behavior of fast electrons in toroidal discharges was per- 
formed. The kinetic equation, describing evolution of fast particle distribution 
was derived and analyzed. Semi - analytical solution of the kinetic equation 
was obtained into suprathermal energy region. External electric field, anoma- 
lous transport, nonuniformity of mean magnetic field and collisions was shown 
to be an important factors, affecting distribution function. The role of am- 
bipolar electric field was established and identified as an essential factor of the 
process of diffusion of fast electrons. The effect of strong influence of density 
profile on diffusion of fast particles was clearly demonstrated. Comparison 
with the experimental data, obtained on ZT -40M device, was carried out. 
The agreements with the results of these experiments is observed. 

1 Introduction 

Plasma confinement in toroidal devices is an important problem of thermonuclear 
fusion research. It is well known that various instabilities leading to discharge tur- 
bulization are readily excited in plasma. Experiments have shown that electron heat 
conductivity is by several orders of magnitude higher than the limit predicted by 
the neoclassical theory |lj . Such a strong energy transfer to a discharge wall is the 
main energy loss channel in toroidal systems. This is the reason why this problem 
is of great interest, see the review 0. 

The theory suggests two possible mechanisms of anomalous transport. One of 
these mechanisms is due to turbulence generated by potential electric field oscil- 
lations caused by drift instabilities 0, ^, [|, |6|. The other mechanism is due to 
turbulence induced by magnetic fluctuations. On exceeding the threshold value 
(determined by overlapping of mean resonance modes), magnetic field fluctuations 



lead to a stochastic wandering of magnetic field lines about the discharge and, as a 
consequence, to a transport of particles moving along these field lines. The qualita- 
tive picture of this phenomenon was first presented in the paper as a stochastic 
description of a diffusion of magnetic field lines. The theory of this mechanism 
of anomalous transport was further on comprehensively studied by many authors 



rPL |TT|, O. According to the theory, the effective diffusion coefficient in the 



former mechanism is inversely proportional to the velocity of particle motion 

D„ ~ <4^ (1) 
Bl v e k ' 

whereas in the latter case the particle diffusion coefficient increases with increasing 
particle velocity 

D ±b ~±-J-L c v e (2) 

The role of each of these mechanisms in the formation of energy flux has not yet 
been finally established. Measurements of the excitation level in tokamaks have 
shown that the first type transport mechanism dominates in the vicinity of discharge 
boundary |]T3], [TJ], [T^|, while inside a discharge (where the level of magnetic 



fluctuations is difficult to measure) the transport may be due to both excitation 
mechanisms. As distinct from the tokamak, in pinch type devices with a reversed 
field, where the magnetic fluctuation level is by two orders of magnitude greater 
than that in tokamaks |T7|], predominant is a transport due to large-scale magnetic 



fluctuations [T^ . 

Along with a thermal particle flux, there also exists a superthermal fast particle 
flux. Moreover, if the particle lifetime in a discharge increases, the number of su- 
perthermal particles may increase appreciably. This effect is obviously of particular 
importance in devices in which plasma is heated in an ohmic way. One of such 
devices is a reversed field pinch (RFP) where, according to experimental data, a 
considerable part of energy is transported by fast particles ]TjJ ^(J, j2~Tf . 

To estimate the energy confinement efficiency in RFP type devices, it is necessary 
to determine the distribution function of fast particles in a discharge. Furthermore, 
fast particles bear information on the processes proceeding inside a discharge, and 
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so the investigation of their distribution may be used to diagnose the state of the 
plasma 0, |22|. 



In the paper J^l , the authors formulated a consistent kinetic theory of anoma- 
lous transport processes in a turbulized plasma and derived the kinetic equation 
describing the averaged particle distribution function in these conditions. A re- 
laxed state of a turbulent plasma and the anomalous transport processes under 
RFP conditions were analyzed in the papers [EHl 23 1. Here we are considering the 
superthermal electron distribution function in RFP. As has already been mentioned 
above, predominant in RFP is an anomalous transport due to magnetic fluctuations. 
Therefore, the influence of this particular mechanism of transport upon the distri- 
bution function of superthermal electrons will be our prime concern in the sequel. 
The anomalous transport produces a strong effect on fast particles because the dif- 
fusion coefficient (|2]) grows linearly with increasing particle velocity. Another factor 
affecting strongly the distribution function of electrons is an electric field applied to 
the plasma. 

It is known that even a weak electric field applied to a plasma induces the 
formation of a tail of runaway electrons in the energy range exceeding the critical 



one, e > e c , where e c = -jf T e , and E c is the critical field |25|, |26fl. In RFP, an applied 
field is high, that is, the ratio E/E c is significantly larger than that in tokamaks, and 
therefore the distribution function distortion due to the action of the field E is much 
stronger here. In addition to the vortex electric field, applied to plasma, a potential 
electric field, a so-called "ambipolar field is generated in the discharge due to the 
difference in the ion and electron coefficients of diffusion. Producing an immediate 
effect upon fast electron diffusion, this field is responsible for the dependence between 
the ejection of fast electrons and the transport of the ion plasma component, which 
have typically been considered independently. On the other hand, as the number of 
fast particles increases, they themselves may start affecting the macroscopic state of 
the plasma, which is for example the case with convective transfer in a rippled field 
[P7|. We shall not consider here the effects due to corrugation, but this role will be 
played by anomalous transport in a turbulent plasma. Thus, an examination of the 
distribution function of electrons is necessary for a correct self-consistent analysis of 
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a macroscopic relaxed state. 

Besides the factors listed above, the distribution function is noticeably affected 
by inhomogeneity of the mean magnetic field and temperature. For example, nonuni- 
formity of electron temperature results in thermal runaway and appearance of hot 



particle tail in the cold plasma region |28 



It should be noted that the influence of anomalous transport in RFP on the dis- 
tribution function of fast electrons was examined in the paper P0 |, but the authors 



proceeded from the model kinetic equation which disregards a number of essential 
factors affecting the distribution function In particular, the effect of Coulomb col- 
lisions, a potential ambipolar field and inhomogeneity of the mean magnetic field 
were neglected. It is therefore necessary to investigate the distribution function of 
superthermal electrons more thoroughly making allowance for the influence of all 
essential factors. This is just the goal of the present paper. In section 2 we derive 
the master kinetic equation with account of the influence of the applied electric 
field, an anomalous diffusion, a potential ambipolar field, collisions and magnetic 
field inhomogeneity. The principal parameters determining the distribution function 
of fast electrons are discussed. Bearing in mind complicity of studying the complete 
problem, at the beginning of section 3 we analyze the distribution function of fast 
electrons in a homogeneous magnetic field. The strong influence of inhomogeneity of 
the electric field applied to plasma and the dependence of the electron diffusion rate 
on the profile of the thermal particle density are resolved. In section 4 we investi- 
gate the influence of inhomogeneity of the mean magnetic field upon the distribution 
function of electrons and point out a substantial deformation of this function. Fi- 
nally, in section 5 we estimate the influence of the indicated effects in specific RFP 
conditions and compare the developed theory with available experimental data. The 
results of experiments are seen to be in agreement with the theory. 

2 The kinetic equation for superthermal electrons 

Let us consider a magnetized plasma with a magnetic field B{r) which has a regu- 
lar B (f) and a fluctuational b components. The amplitude of fluctuations will be 
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assumed small as compared with the mean field Bo, \b\ <C |-Bo|- The basic quanti- 
ties characterizing the fluctuations (the correlation length and the correlation time) 
will be thought of as large as compared to the Larmor radius of particles and their 
inverse gyrofrequencies as in |22j. Within such a statement of the problem and dis- 
regarding toroidality effects, the authors derived the general kinetic equation for the 
distribution function of particles f(r, u, fi) averaged in the ensemble of fluctuations, 
which, by virtue of cylindrical symmetry about the angle 9 and the direction z along 



the cylinder axis depends only on the radius r [23 



% + — eM = St{f) + I{f) (3) 
at m e ou 

where the collision integral of particles with fluctuations is given by 

'</) = w) + - lw iw 

K _ H Fdf Hdfrj^Ea _ hdBU 
u BW + \u\~du i rn e B ?WBJ 

Here the ambipolar electric field eE a = ^ + 37^) {1 ~ 2<5 m } due to the differ- 
ence between the diffusion rates of electrons and ions, the constant S m = ^Jm e /mi <^ 
1 gives correction related to ion diffusion (assuming that transport of ions is defined 
by magnetic fluctuations too), n(r) is the particle number density, F = J °° dL(b r b' r ) 
-is the correlation function of fluctuations b r , the integration over L goes along the 
trajectory of particle motion, b' r = b r (r(L),9'(L),z(L)), u is the particle velocity 
along the magnetic field line, ft = u\/B is the adiabatic invariant of particle motion, 
E e is the external field, E e = E e ■ h, h = B Q /B Q , St(f)- is the Coulomb collision 
integral of particles. 

In plasma heating devices, the number of fast particles Nf is always small as 
compared to the concentration iV of main particles. That is why the kinetic equation 
(|3[) can be linearized in the small parameter Nf/N 1. Equation (|3]) is convenient 
to write in a spherical coordinate system in the velocity space e, /x, where e = u 2 +u 2 ± 



is the total energy and \x = u/y u 2 + u\ is the cosine of pitch angle (below for brevity 
we'll use simply pitch-angle without cosine). Linearizing we obtain 



% + ^(r)5 l{ 2,e% + (1 - , 2 ) d M = 4T(r)U + 2% + 



2E a rdK , l~lJ?dK x dB 1 l-^idK x m 

Here E e (r) is the profile of the external longitudinal electric field normalized to 
the critical electron runaway field £7 C E3], £? c = 47re An ^ j£ e ^ is the effective ion 

e 

charge, A is Coulomb logarithm, Si = E e o/E c is the dimensionless parameter char- 
acterizing the magnitude of the longitudinal field E £ q relative to the critical field E c , 

S2 = v a /i>Q is a dimensionless parameter characterizing the particle-fluctuation col- 

F I T 

lision frequency v a = y ^f- as compared with the Coulomb collision frequency 

of electrons z/ = ^1% n i) 2 ? a * s the characteristic system dimension, T(r) - is the 
mj TJ 

profile of the electron temperature normalized to the temperature at the center T e . 
Furthermore, dimensionless quantities r = r/a, F = F/F max , B = B/Bq, E a = 
E a ea/T e , e = e/(T e /m e ), r = v§t are introduced, the sign ~ marks the correspond- 
ing dimensional quantities, S\ and S2 are small parameters of our problem. In what 
follows we shall consider steady-state solutions of equation (|) to determine the 
established distribution function of superthermal particles. We shall assume the 
distribution function of the main particles in plasma to be stationary and equilib- 
rium. In such a statement, it will be a source of superthermal plasma particles. As 
boundary conditions for the distribution function / it is natural to require that / 
be regular as r — > and that all the particles die on the boundary for r = 1, that is, 

^| r=0 = o /U = o 

The investigation of the complete problem is difficult because of a simultaneous 
action of such factors as an applied electric field E e , an anomalous diffusion and 
inhomogeneity of the mean magnetic field. We shall therefore begin with examining 
a joint effect of the external electric field E e and the fluctuations assuming the 
magnetic field gradient to be small 

dB a 

~dr~B 1 (6) 
and then proceed to the case of an inhomogeneous magnetic field B(r). 
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3 Distribution function of superthermal electrons 
in a homogeneous magnetic field. 

Provided that the condition (^) is satisfied, equation (f|) has the form 
E e (r)h{2m% + (1 - = 4T(r)0 + 2% + 

Ba = (f ^ + ^){1 - «, S m « 1 

Equation @ describes the established distribution function of electrons in the pres- 
ence of the field E e and plasma turbulence. The method of solving equation (H) 
depends on the energy range within which we seek the solution. Therefore we shall 
first examine the energy range immediately adjoining the equilibrium region, which 
is henceforth referred to as a polynomial region of solution 

1 < e < ^ 1/2 , 8 X < 1 (7) 



We shall not consider here the thermal runaway of particles [^3, According 
to the results of |29[ one may estimates the critical value of energy y^, y = e5/ 
when the temperature profile will relax to homogeneous one, yu = 

ll c 

where a - is a scale of the system, L\\ c - is a correlation length of fluctuations along 
magnetic field line. As usually -j^- <C 1 and we assume that, 

Li\\c 

Vk < 1 (8) 

And therefore, in accordance with @ in equation(|6]) we have put T e = const and 
for the sake of simplicity Tj = const, because ion temperature profile contributes 
only to a correction term in ambipolar field. So, as is seen inhomogeneity of electron 
temperature contributes only to ambipolar field. 
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3.1 The behavior of solution in the polynomial region 

To begin with, we consider the case of a constant field E e (r) = const. We imply 

I/T-l 

at first to clear up the effects of ambipolar field, that *C S m and neglect 

contribution from temperature gradient in ambipolar field. In equation (^j) we pass 
over to a new variable y = eo/ and seek the solution as a series of eigenfunctions 
of the Sturm-Liouville problem 



dxi 
dr 

dxi 
dr 



that is, 



r=o — 

r=l = 
CV2 



/ = Cn(r)exp(-v]/ / ( 5i /2 )^ X/ (r)^(z/,/i) (10) 



As the boundary condition for y — > we require the condition of sewing with the 
equilibrium distribution function 

fo = Cn(r) exp(-y/5l /2 ) 



n(r)| r= i = 

-1 /r) 

Substituting (|10|) into (|6|) and collecting terms with the same power 5/ in zero 
approximation, we obtain ty = y/2. In the next order of perturbation theory, 
neglecting correction of the order of S m <C 1 in the expression for the ambipolar 
field (H), we obtain a system of equations for the functions Ri(y,fi) 

2 f -T i |« 1 -^f>-*l^' + 2 "^ < u > 

with boundary conditions as y — > 

i? (0,/i) = l, i2/(0,/x) = Z^O 

> 

where (3 = 62/ Si. Since zero eigenvalue of the problem (H) ?7o = 0, for the function 
Ro(y,fi) we obtain an equation containing no contribution of anomalous diffusion, 

8 



and it will have a solution similar to the one obtained in the paper p5| , where the 
distribution function was distorted only by the electric field E e and had a directed 
character for fi ~ 1 (Fig. 1, the dashed line). The solution of the system (|TTD for 
/ 7^ 0, which satisfies the boundary conditions (because for I ^ all the eigenvalues 

-1/2 

r/i > are nonnegative), will be Ri(y,fi) = 0. So, up to terms of the order of 8 X 

1/2 

(terms of the order of 5/ 5 m <C 1 and 5i <C 1 are neglected) the initial equilibrium 
distribution function f will not be distorted by anomalous transport. This means, 
as is readily seen, that the contribution from the anomalous diffusion is completely 
compensated by the ambipolar field E a in exactly the same way as in the case of 
diffusion of thermal plasma particles. As is well known, electrons and ions diffuse 
together as a single whole with a doubled ion diffusion coefficient, which in our case is 
5 m times smaller than the electron one (provided that the predominant ion diffusion 
mechanism is also an anomalous transport caused by magnetic fluctuations). 

To determine the distortion of the distribution function by an anomalous trans- 
port, we have to examine the solution of equation @ with allowance for terms of 

1 /2 

the order of 8-y 8 m and Si , which under certain conditions (to be discussed below) 
become appreciable. Now we are in a position to consider the case where the param- 

1 12 

eter 8 m satisfies the condition 8 m ^> Oj/ . This means that we shall take into account 
the contribution from the anomalous diffusion, proportional to the anomalous ion 
diffusion coefficient. 

With allowance for corrections of the order of 8 m , the solution of equation (|6|) 
is convenient to seek as before in the form of an eigenfunction series of the Sturm- 
Liouville problem 

^{rF^} + X l9l = (12) 



dm 

dr 



r=0 — 

gi\ r =i = o 

1/2n 



f = C exp^o/^ ) &ft(r)e,(y,/i) (13) 
n{r) = EiAgi 
with the boundary condition as y — ► 0: 

e,(o,//) = A, (14) 



In order that we might pass over to a system of ordinary differential equations, 
it is also convenient to expand the functions Oi(y,n) in a power series of Legendre 
polynomials P m (ji): 



Substituting flT3| ) into (|B|) with account of flI5p, we ultimately arrive at a system 
of linking equations for the functions F™(y), = y/2: 

rlF m Z t * 

2 ^T = ^ m(m + ^ + y^HmnP + (16) 

pra-X ( m+ l)jTm+l 

Zy {'2m-l + 2m + 3 > 

with boundary conditions 

if (0) = 4 
F; m (0) = m ^ 



where 



= Io r 9i{f-§f{ rIl k} ~ E a U k }dr/ Jorgfdr 
U k = F{^-E a g k } 
'j mn = (2m+l) j\P m P n \jj\dn, p = S 1 /S 2 



The system of equations (fig ) is solved numerically by cutting off the chain of 
equations on the term L in the expansion in gi and on the term M in the expansion in 
P m , so that a doubling of L and M changes the solution by less than 10 %. Clearly, 
strong distribution function distortions by an anomalous transport will take place 
only when the parameter A = 2f3X m 5 m is of the order of unity A ~ 1, where A m is 
the maximum eigenvalue corresponding to the eigenfunctions gi(r) which contribute 
to the expansion (|i4|). The smallest distortions may be expected in the case n(r) = 
9o(r). 

We are examining the behavior of the system (|T6|) in the model case F = const, 
when the eigenfunctions of the problem ( |T2"D are the Bessel functions gi(r) = Jo (£/?")• 
We begin with examining the case n(r) = go{r). The solution of the system fll6|) 
depending on the pitch angle fi and the energy y is presented in Figs. 1,2. Figure 1 
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shows solutions obtained for various values of the parameter (3. For small (3, when 
A < 1 and the effect of the electric field E e dominates, the solution is close to that 
obtained in the paper [25| (the dashed line in Fig. 1). As (3 increases, the influence 
of the anomalous transport (when A > 1) becomes predominant. In this case the 
solution becomes almost symmetric in /j. As should be expected, it is concentrated 
in the vicinity of /x = , where the diffusion coefficient (0) vanishes and falls 
symmetrically for fi ±1, where the diffusion coefficient is maximal. Depending on 
the energy (Fig. 2), the distribution function falls exponentially ln(/// ) ~ —A y 2 , 
where A = 2S m (3\o (Ao is a nonzero eigenvalue of the problem (|T2|)). The function 
/, which is shown in fig. 2, was averaged over //. 

It should be emphasized that the initial distribution over discharge n(r) will not 
be deformed considerably. 

We shall now see what will happen if as the initial profile n(r) we shall choose 
an arbitrary function n(r) such that the expansion (|T1| ) will involve (and with a 
substantial contribution) harmonics with I ^ 0, as is shown in Fig. 3. The harmonics 
with eigenvalues \i ^> A should be expected to damp faster than the zero component 
does already for y < 1. As a result, the distribution over discharge will rather 
rapidly relax to go{r), which will immediately lead to balance violation between 
the ambipolar field and anomalous diffusion. In this case, an effective increase 
of the diffusion rate may be expected. Indeed, Fig. 3 shows the distribution over 
discharge radius for y — 1, obtained in the solution of the system (|T6| ) with the 
profile n(r) = n(r), is almost coincident with go(r) (dashed line in the same figure). 
On the other hand, Fig. 2 presents the dependencies obtained for one and the same 
value of (3 for n(r) = g(r) and the profile n(r) = n(r) depicted in Fig. 3. It is seen 
that in the latter case the decrement is considerably larger, which testifies to an 
effective increase of the diffusion rate. 

The result was obtained make allowance us to conclude, that the rate and charac- 
ter of diffusion (and as a consequence lost of energy due to diffusion of fast particles) 
essentially depends on profile of mean particle density. The minimum of lost will be 
as n (r) = g (r). 

The density profile in a device n(r) is defined by many factors: anomalous trans- 
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port, neutral particles flow from the camera wall and external sources, accelerated 
ions injection, convective transport processes. So by means of injection, for example, 
one may to drive profile n(r) and respectively process of fast particle transport. 

The result obtained suggests that the compensation of anomalous transport by an 
ambipolar field is a consequence of equilibrium in the leading term of the distribution 
function of electrons, that is ^o{y) = y/2, and the absence of distortion of the initial 
electron distribution n(r). As shown above, the strict balance is violated by the 
anomalous transport itself when n(r) ^ (?o( r )- 111 the presence of an electric field 
such a balance will also be violated for y » 1 since according to |2I| the electric 
field E e induces strong deviations of electron distribution from an equilibrium one. 
We may point out another mechanism of balance violation already for < y < 1, 
namely, inhomogeneity of the applied field E e (r) which in this case plays the role of 
a fast particle source nonuniform in space (note that it is exactly the case realized in 
RFP). We shall consider the solution of equation (|6]) in the polynomial region using 
the methods presented above. As before, we put F = const. Figure 2 presents the 
dependence of the distribution function ln(/// ) on the energy y for one and the 
same value of the parameter f3 and n(r) = go(r) in two cases: 

a) for E = const and b) for E e (r) = Jo(n r), k = 3.0, (such a profile is close 
to the one observed in RFP discharges). It can be readily seen that in case b) the 
distribution function fall with energy is much higher than in case a), which shows an 
acceleration of the diffusion process. It is noteworthy that the electron temperature 
inhomogeneity (which we do not consider here) will obviously play a role similar to 
that played by the inhomogeneity of the external field E e (r) and will also lead to a 
violation of strict balance, see Fig. 2. 

We have assumed above that 5 m ^> <V , and the correction of the order of 

1 /2 

5 m 5i has led to strong distribution function distortion when A > 1. Clearly, 

1 /2 

in the converse case, that is when 5 m <C <5 X , there will proceed an analogous 
process which is not distinct qualitatively from the one investigated above with the 
only difference that in this case the role of the parameter A will be played by the 

1 /2 

parameter A x = 2/3A m <5 1 . As a result, as before for Ai ~ 1 the distribution 
function will relax rapidly to the profile (?o( r )- 
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Thus we have completely investigated the behavior of the solution of equation (H) 
in the polynomial region. We have obtained that in a special case where the profile 
n(r) is chosen in the form n(r) = go(r) and the external field E e is homogeneous, 
the ambipolar field E a completely damps the anomalous transport of fast electrons 
with an accuracy of 6 m <C 1, so that the effective diffusion is determined by the 
parameter A or Ax which is substantially smaller than f3X , A <C f3X . On the other 
hand, such a strict balance in higher terms of expansion may be violated provided 
that n(r) ^ go{ r ) or the external field E e and temperature T(r) are inhomogeneous. 
We note that the latter always takes place in RFP. An essential result is here the 
fact that for A > 1 (or Ai > 1) the distribution function of fast particles over 
discharge relaxes with increasing energy to the universal profile (?o( r ) independent 
of the initial distribution. This fact will simplify appreciably our analysis in the 
remaining part of this section, where we consider the energy range y » 1. 

3.2 The behavior of the solution in the exponential region. 

We shall now consider the domain of solution for high energy values e ^> <5 -1//2 , 
where 6 = 61,62- The distortions of the distribution function in this domain are 



known to be of exponential character [26]. Therefore, a polynomial expansion is not 



effective here. As we have seen above, practically for e 3> 6 l l 2 , the distribution 
function of fast electrons over discharge is coincident with the zero eigenfunction 



go(r) of the problem (12). Therefore, in this energy range, in the expansion of the 



distribution function it is natural to make allowance only for terms containing go(r): 

1 

Let us consider the case E e — 0. We shall pass over to a new variable z = e<5 2 
and represent \I/ as: 

^ = V /6 2 + ^i/5 2 1/2 + * 2 + . . . (18) 

_____ ^2 

Substituting (|17D and (|18|) into (|) and collecting terms with the same powers 6 2 , 



we obtain the system of equations 

9*0 



3, ° (19 > 
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(20) 



4 ^ f _^ + ^ {l _^_ :M{Xa + B ^ } ^ (21) 

B = -4 / rgl(r)F(r)E 2 a (r)dr/ [ rgl{r)dr 
Jo Jo 

From (|19D it follows that — ^0(2), an d ([^) holds identically. From the condition 
of the absence of a jump of the derivative for /x = there follows a natural 
condition (in view of symmetry under a substitution of /x for — /x in (^) for i? e = 0) 

mi \ -0 

Taking into consideration this condition, as well as the fact that \l/ = ^?o{ z ), from 
equation ( 21]) for /x = we obtain 



Substituting the expression found for \I>o (-^) back into (plj) , we obtain the equation 
for determining ^i(z, /x). 



|/x|A* 
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<9/X \| Z e// (1 - /X 2 ) 

where Aq = Aq + -B/4 is an effective eigenvalue with account of the influence of 



the ambipolar field E a . From ( 22"D we see however that the asymptotical expansion 



obtained is violated in the vicinity of /x = because the second derivative ^ % 

contains a singularity. Indeed, ^% ^ oo /x 0, which indicates of the presence 
of a boundary layer near /x = 0. To obtain a correct expansion, it is necessary to 
investigate the behavior of the solution in the vicinity of /x = 0. To this end, the 
small term with a second derivative ^ % b\[ 2 in equation (|21|) should be retained. 
Omitting terms of the order of /x 2 as /x — > in ([H]) and making a substitution 

1 /2 

= — ln(0)5 2 , we come to the Airy equation 

^ = (23) 

14 



^)-«^-^) 

Equation ( ^3|) has two linearly independent fundamental solutions, one of which 
grows exponentially as £ — > oo and may be discarded for being limited. The other 
solution = CqAi(£) is shown in Fig. 4, £ is the point where the function O and 
has an extremum, = 0. Then from the condition | M=0 = we obtain 

an additional relation that allows us to determine the function ^fo(z) 

4(^) 2 - 2^° = Uzf f X^\z)^z^ (24) 

In the general case, to find ^o(z) it is necessary to solve the nonlinear equation 
(p4|) ) but for z <C S^ 1 expanding equation (^4j) in the small parameter §2 we obtain 

*oW = f + ||&|^AS 2/s *i /S ^ s (25) 
A* = A + B/4 

Taking into account (|2"5]), we may determine the angular dependence of the distri- 
bution function: 



^{A-^^l^l-ieol^A-^) 2 / 3 ^ 3 ^ 3 } ^ 
\ r, 2. dfi + *i (z) (26) 



where ^i(z) is an unknown function. As fi — > 0, it is necessary to sew fl26|) with 
the solution of equation fl23|) to find the constant C . It is readily seen that as 
$2 ~ *■ 0, (|25"D transforms into ^o(z) = z/2 and ( p6|) into fl22|) . As z — >■ 0, the solution 
found here must pass over to the solution obtained in the polynomial region. Figure 5 
shows dependencies of the distribution function on the pitch angle in the polynomial 
region for y > 1 and the solution obtained in the exponential region when z — > 0. 
One can see good agreement between the two solutions. The solution obtained by 
us shows that the leading term of the asymptotical expansion for 62 — > is close 
to equilibrium, and deviations from it occur only for z ~ 5^ ■ This result is in 
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close agreement with the results of the direct numerical simulations pl| carried 



out for TOKAMAK conditions in the absence of the applied electric field E e and 



experiments [ 32 1 , [33] . It should be noted that the role of an ambipolar field in 
this range of energy values comes down to renormalization of the eigenvalue Ao 
which determines the anomalous transport, namely, the effective eigenvalue with 
allowance for the influence of the ambipolar field E a will be Aq = Ao + B/4 , Aq < Ao 
since B < 0. That is, as expected, the ambipolar field E a damps the anomalous 
transport of fast electrons. It is clear that, this effect is important when the effective 
temperature of fast electrons is of the order of equlibrium temperature at the center 
of the discharge. 

The field E a also affects the behavior of the distribution function in the far region 
of energy values as z — > oo. So, when B ^ one can readily obtain an asymptotical 
expression for ^q{z) as z —* oo: 



At the same time when B = and z — > oo, from (|24|) we have 

3 

14 



* (z) - ^ 7/ % /3 4 6 A /e 



that is, in the absence of the ambipolar field E a the distribution function falls 
stronger with energy. 

We shall now consider the case with a nonzero external field E e . As in the 
previous case in the expansion (|T7|), we shall take into account only the contribution 
from the terms go{r). In equation (P) we pass over to a new variable z = e5\ and 
represent the index of the exponential in (|PTD in the form 

m = + m^s 1 / 2 + * 2 + • • ■ (27) 

Substituting (|27| ) into @ and collecting terms with the same powers S^ 2 , we obtain 
a chain of connected equations for the functions ^O; ^i, ^2, ■ ■ ■ '■ 



d/j, 

9^o 
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(28) 



(29) 



4(^) 2 - 2^ + Z -f{l - - zW n*) + 2^ = (30) 



<9* n <9*i <9^1 n 9*1 

^anr^r^ 1 -^' 1 (31) 

B = -A I rgl(r)F(r)E 2 a (r)dr / [ rgl(r)dr 
Jo Jo 

E = j rgl{r)E e (r)dr I \ rg%(r)dr 
Jo Jo 

P = 6 2 /6 1 ; \*( Z ) = \ + B(^f 

oz 

From the first equation (f28| ) it follows that *o = *q(<z). Given this, equation 
(|29| ) holds automatically. Then, in the absence of a singularity for \i — 1, from (|3~0D 
we obtain the equation with the help of which we may determine *o : 

4(-^— ) -2— - z(3\ (z) +2zE— = (32) 



This implies 



(1-zE) + J(1-zE) 2 + z[3\ (4-PzB) 



dz 4 - 2/35 

Integrating (133) we arrive at 



(33) 



5^ 5^~ 



(8a 2 + a\Bj3) ln(ai + 2a 2 £ + 2^fa~2\j 1 + + a 2 -2 

(45 - 5/3) ln(gi) 
5^ 



5i = -4 ai 5 3 /3 3 - 25 4 /3 4 - 8a 2 B 3 (3 3 z - ai B 4 f3 A z - 25 3 /3 3 (45 - 5/3) x 



(34) 



x a/1 + a\Z + a 2 2; 2 
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Oi = 4(3X - 2E- a 2 = E 2 - (3 2 X B 
Taking into consideration Q3"3|), we obtain from 



/, /3(A + 5(^) 2 )(|/i|-l)-2^(/i-l) 

9 l {n,z) = -z . + (35) 

A ^ z eff (l-n ) 

It should be noted that holds for 

f3{\ + 5/4} < 5 (36) 

In this case, the solution will not have a singularity for \x = 1. Substituting (|35| ) into 
, provided that there is no singularity for fi = 1, the unknown function ^1(2;) 
may be determined: 



= ^("Vv^ + 1) f L| m ^—dz (37) 

Jo 4—^ + z(3B— ± + Ez-1 
oz oz 

The expressions (|34]), fl35|) and (|57|) determine in the exponential approximation 
the dependence of the distribution function in the exponential region on the energy 
and the pitch angle /1. As z — > 0, the obtained solution must coincide with the 
solution in the polynomial region for y > 1. Figure 6 shows both the solutions 
which are seen to be almost coincident. 

To establish the asymptotical behavior of the distribution function as z — > 00, 
we represent the index of the exponential in fllTD in the form 

M> = \fr QZ + ^r x + . . . 
^0 = <Po/S! + <Pi/5l /2 + <p 2 + ... (38) 

z = edi 

Substituting ( [38f) into @ with allowance for ([FTP , keeping terms of the order of 2 

1/2 

and collecting terms with the same powers <V , we obtain the system of equations 
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for ip , ip 1} . . .: 



(39) 



{^lfZ eff {\ - - H/3{A + Btf] + 2^ E = (41) 
2Z eff (l - /i 2 )^^f - WBcp oVl + 2^E Vl + (1 - /x 2 )^£ = (42) 

where (3 = 82 /Si, , and B and E are similar to fl30D , whence ipo = const. Then in 
the absence of singularity for fi = 1 we have from fl4"T|): 



^0 = ^ (43) 



Knowing ip from fl4"I|), we obtain 



eta 



n(4E/f3B)(JE 2 - /3 2 X B - E) 



<fi = const jj, > 

nUF.Il 

<Px = ~ Ii 



fi<0 (44) 



^4E/(3B)(^E 2 -(3 2 X B-E) ^ 
\ t^t, ~Ts dp + ¥1 // < (45) 



where (pi is an unknown constant. Then in the absence of singularity we obtain 
from (H) for /i = 1 that ^ = 0. 

The expressions ( [4"3"D and ([EJ) determine the asymptotical behavior of the dis- 
tribution function for E e (r) 7^ as z — > 00. It can be readily seen that in the limit 
z — > 00 (0) goes over to ( f43|) and ( p5|) into (f4"5|). So, the obtained solution 



( p5|) in the exponential approximation describes the particle distribution in the en- 

1 /2 

tire range of energy values for z > 5 l . It should be noted that the asymptotical 
expression for the distribution function for z — > 00 (f£3|), fl4"5] ) holds also under the 
condition inverse of (^), but obtaining this asymptotes is apparently not described 
by the simple expansion (|i~Tp. The solution obtained in this section describes com- 
pletely the distribution function of fast particles in the presence of an external field 
E e and anomalous diffusion in a homogeneous magnetic field. We have investigated 
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here a rather general case in the model assumption F = const. In section 5 we 
shall analyze the distributions of fast electrons in specific experimental conditions, 
and now we proceed to the question of the influence of inhomogeneity of a mean 
magnetic field. 

4 Distribution function of fast particles in the 
presence of a finite magnetic field gradient. 

We are now ready to consider the influence produced by an inhomogeneous magnetic 
field upon the distribution function of fast electrons. From the point of view of 
physics, the occurrence of terms proportional to the magnetic field gradient in the 
kinetic equation (|3]) is associated with keeping the adiabatic invariant u\/B{r) = 
const upon a diffusion particle motion along a discharge. In this case, particle 
diffusion is responsible for particle energy redistribution from the transverse degree 
of freedom to the longitudinal one and vice versa, depending on the sign of the 
derivative dB/dr. Taking into consideration the fact that on the average particles 
move from the center to periphery of a discharge, in the case dB/dr < part of 
the transverse particle energy is transferred into the longitudinal and conversely in 
the case dB/dr > 0. Such dynamics must in turn affect the diffusion process itself. 
Indeed, in case the amount of particles with high longitudinal energy increases (when 
dB/dr < 0) or on the contrary decreases (when dB/dr > 0) because the anomalous 
diffusion coefficient (0) is proportional to the longitudinal velocity, the diffusion rate 
will respectively either increase or on the contrary decrease. 

Let us consider the stationary solution of equation (|3]) in the range 1 < e < <5 -1 / 2 , 
where 5 = 5i, 5 2 . Not to complicate the picture, we put the external field E e {r) = 0. 
In order to obtain the solution, we shall use the technique we applied in the solution 
of equation ([]) in this range, that is, the polynomial expansion (|i"5|). The function 
F(r) will put F = const, and n(r) = go(r) of the problem ( fl"2|) . Figure 7 shows the 
dependencies of the distribution of fast electrons obtained in the solution of equation 
@ depending on the electron energy y in two cases: dB/dr > and dB/dr < 0. 
The figure also presents for comparison the result for dB/dr = 0. As has been 
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expected, the greatest decrement of ln(///o) appears in the case dB/dr < 0, and 
the smaller in the case dB/dr > 0. It would also be of interest to trace the behavior 
of the solution obtained depending on the longitudinal u| and the transverse u\ 
energies. In Fig. 8 we can see the dependencies of ln(/// ) on the transverse and 
the longitudinal particle energy in the cases dB/dr < 0, dB/dr > and dB/dr = 0. 
We see that both for dB/dr = and dB/dr > there exists anisotropy in particle 
distribution over transverse and longitudinal energy, so that T± > Ty. However 
for dB/dr < the anisotropy changes sign and then T± < T», which is naturally 
connected with a substantial energy redistribution from the transverse degree of 
freedom into the longitudinal one. This effect is most clearly pronounced in the 
dependence of the distribution function on the pitch-angle /x, Fig. 9, which differs 
substantially from that obtained in the case dB/dr = 0, Fig.l. The fall in the center 
(for /i ~ 0) in Fig. 9 and the appearance of maxima near fi ~ 1 is just due to particle 
efflux from the region /i ~ towards the region /i ~ 1. 

5 Behavior of the distribution function of fast 
electrons in RFP. 

In the preceding sections of the paper we have analyzed the influence of various 
factors upon the distribution function of fast particles. It should be noted here 
that just in RFP discharges the contribution from each of them will be substantial 
because all of them are present in this type of devices. An RFP discharge possesses 
a force-free magnetic configuration with components of the mean magnetic field (the 
toroidal B z (r) and poloidal B$(r)) which are of the same order of magnitude and are 
determined by one and the same parameter = 27raI/Q, where / is total current, a 
is small radius of the torus, and $ is the total toroidal magnetic field flux. Magnetic 
field components are well enough described by Bessel functions of the form 

B z {r) = B J (2Qr) 

(46) 
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B e (r) = B .h(2Qr) 

see Fig. 10, [34]. So, the scale of the magnetic field gradient in RFP will be of the 
order of the small radius a of the torus : 

dB a 

~ 1 

dr B 

Since plasma in RFP is heated in an ohmic way, present in the discharge is 
a sufficiently strong longitudinal electric field E e . Since an applied electric field 
E = Ee z is toroidal, its projection E e (r) onto a magnetic field line is described by 
the relation 

Ee(r) = E^f\ (47) 
B{r) 

From fl4"7|) it becomes clear that the applied field is strongly inhomogeneous about 
the space. 

Direct measurements of distribution of superthermal electrons in steady - state 

of RFP were carried out on ZT-40M devices |fi9| , |35| , |36| , j37| . For comparison with 

the theory, we shall consider the experimental data obtained in P5fl . Performed in 

the experiment were direct measurements of energy distribution of fast electrons 

in the near-boundary region of discharge with a temperature T w ~ 20eV near the 

boundary and To ~ 220eV in the center of discharge. The characteristic value of 

the parameter = 1.4. Observations were carried out for particles of energy up to 

1.5keV. Against the background of cold plasma with temperature T w ~ 20eV one 

could see a tail of energetic electrons which moving in the direction of the magnetic 

field and have Maxwell distribution with characteristic temperature 7\ ~ 530eV and 

a flux of particles that moved backward and had a temperature T 2 ~ 330eK. The 

anisotropy in the distribution of particles in the longitudinal and transverse energies 

was observed, the temperature being T± < Tji. Under conditions of the experiment 

described above, the parameters Si, S 2 of our problem were 5i = 0.21, 5 2 = 0.036 for 

To = 220 eV, n = 2,4 • 10 13 cm~ 3 , E e = lOV/m and the amplitude of fluctuations 

\b/B\ ~ 1%, ~ 4- 10- 4 . 

BqCl 

The solution we obtained in section three (formulae (P4"|), (P7|)) gives qual- 
itative agreement with experiment: exponential distribution of particles on energy 
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and anisotropy in particles distribution over forward and backward direction along 
magnetic field lines. 

An important comparison with experimental data is given by modeling su- 
perthermal electron current I eea , collected by the electrostatic energy analyzer 
(EEA), as function of applied retarding potential V. The current is related to the 
electron distribution function f(r, //, e) by expression: 



oo roo 



Ieea(r,V) oc J f (r, fi, e)v±dv ± V\\ dv\\ (4f 



where 



v = yj2eV/m e 

Figure 11 shows the dependencies of current I eea calculated by formulae 



(p5|), (|37|) , (|4q) for the given values of the parameters 8± and 82 for particles flying 
along the direction of the magnetic field and backward. The profile of the correlation 
function F(r) shown in Fig. 12 was borrowed from ref. [BH]. The graph of the zero 



eigenfunction of the problem ([12]) go(r) for a given F(r) is shown in Fig.13, where 
Ao = 1.2. The profile n(r) was chosen in the form n(r) = go(r) and T(r) = 1 — r 2 . 
The calculated current I eea was normalized on experimental value. The absolute 
value of I eea seems to be rather uncertain because of a lack of information on the 
correlation function near the boundary. 

Figure 11 testifies to not only qualitative agreement between the theory and ex- 
periment (Maxwell distribution and the anisotropy observed in the distribution of 
particles flying in the direction of the magnetic field and backward), but to quan- 
titative one as well. So, the characteristic temperature T{ eor for particles moving 
along the magnetic field makes up T\ eor = 526 eV. But for particles moving in the 
backward direction it was T^ 01 " = 183 eV, i.e. less than experimental one. Backward 
flowing component of fast electrons is seen to be a much smaller than forward one 



(less than 10 % ). This value is also close to experimental observations [35]. 

The behavior of distribution of fast electrons on longitudinal and transverse 
energy is shown on fig. 14, at the parameters mentioned above. In this case for 
particles moving along the magnetic field line T± < T» and 71 1 = 513eV, Tj_ = 385eV. 
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6 Summary 



In conclusion we shall present the main results obtained in this paper. In the 
entire energy range above thermal energies, we have constructed the solution of 
the kinetic equation describing the distribution function of fast particles in the 
presence of an external electric field, collisions, and anomalous diffusion due to 
magnetic fluctuations. The ambipolar electric field resulting from the difference in 
the ion and electron diffusion rates is shown to play essential role in the process of 
diffusion and to be responsible for the fact that the character of the diffusion depends 
strongly on the density profile of the main plasma particles and on the profile of 
an external electric field applied to the plasma. The influence of inhomogeneity of 
the magnetic field of a discharge upon the distribution function of fast particles is 
investigated. The analysis performed suggests that the magnetic field inhomogeneity 
is an important factor affecting the distribution function of fast particles. In the 
polynomial region of energy values, which is an immediate neighbor of the thermal 
region, the semi - analytical solution is constructed with allowance for all the factors 
mentioned above. On the basis of comparison with experimental data we have 
shown that the theory describes the basic details of the distribution function of 
fast particles. We have referred our analysis mainly to the case of RFP discharges, 
but the results obtained in the absence of strong electric field seems to be closely 
connected with anomalous process in TOKAMAK discharges. 
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Captions 



Fig.l The dependencies of fast particles distribution function on pitch angle /i into 
polynomial area of solution under various parameters (5, y = 2, r = 0.5. Dashed 
line - solution, obtained in []25J under (3 = 0. 

Fig.2 The dependencies of fast particles distribution function ln(/// ) averaged 
over /i on energy y into polynomial solution area. 1) - n(r) = go{r), T(r) = const, 
E e (r) = const 2) - in case of nonuniform electric field n(r) = go{r), T(r) = const 3) 
- n(r) = n(r), plot h(r) is shown on figure 3, T(r) = const, E e (r) = const and 4) 
E e {r) = const, n(r) = go(r), T(r) = 1 — r 2 . (3 — 0.3, r = 0.5. 

Fig. 3 The distribution over discharge of initial density profile n(r), distribution 
function f(r) at y = 1, obtained in solving the system ( |I~6D and eigenfunction of the 
problem (|12|) (dashed line). 

Fig. 4 The solution of equation ( p3|) function of £. 

Fig. 5 The fast electrons distribution function f / fo as a function of /x in the 
absence of applied electric field E e . Solid line - the polynomial solution at y > 1, 

1/2 

dashed line - the exponential solution at z = yd 2 ^ 1. 

Fig. 6 The fast electrons distribution function f / fo as a function of u with the 
applied electric field E e present. Solid line - the polynomial solution at y > 1, dashed 

1/2 

line - the exponential solution at z = y5 2 1. 

Fig. 7 The dependencies of fast electrons distribution function, averaged over 
/i, ln(///o) on energy y with finite gradient of mean magnetic field present. 1) 
dB/dr > 0, 2) dB/dr = 0, 3) dB/dr < 0. r = 0.95 

Fig. 8 The dependencies of fast electrons distribution function ln(/// ) on lon- 
gitudinal (solid line) and transverse (dashed line) energy, in the present of finite 
gradient of mean magnetic field. 1) dB/dr > 0, 2) dB/dr = 0, 3) dB/dr < 0. 
r = 0.95 

Fig. 9 The dependencies of fast electrons distribution f / f on u under various 
positions over radius. dB/dr < 0. 

Fig. 10 Radial distribution of magnetic field components across the minor radius. 
Fig. 11 The comparison of theoretical versus experimental current I eea , collected 
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by electron energy analyzer as a function of retarding potential V for the case of 
ZT-40M. Solid line - for particles moving along the magnetic field line, dashed line 
- for particles moving backward. Experimental data are shown by labels. 
Fig. 12 Correlation function of fluctuation in RFP. 

Fig. 13 Eigenfunction of the problem (|I2|) for the correlation function, depicted 
on figure 12. 

Fig. 14 The dependencies of fast electrons distribution function ln(/) on longitu- 
dinal (solid line) and transverse (dashed line) energy for particles moving along the 
magnetic field line. /3 = 0.17, <J a = 0.21, 5 2 = 0.036. 



29 



